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Abstract. High resolution data on low energy phonon branches (acoustic and soft optic) along the three 
principal symmetry axes in pure KTaOs were obtained by cold neutron inelastic scattering between 10 and 
80 K. Additional off-principal axis measurements were performed to characterize the dispersion anisotropy 
(away from the <100> and <110> axes). The parameters of the phenomenological model proposed by 
Vaks pS] are refined in order to successfully describe the experimental low temperature (10 < T < 100 
K) dispersion curves, over an appreciable reciprocal space volume around the zone center (|q| < 0.25 rlu). 
The refined model, which involves only 4 temperature-independent adjustable parameters, is intended to 
serve as a basis for quantitative computations of multiphonon processes. 

PACS. 77.84. Dy Niobates, titanates, tantalates, PZT ceramics, etc. - 63.20.-e Phonons in crystal lattices 
- 62.20.Dj Phonon states and bands, normal modes, and phonon dispersion - 78.70.Nx Neutron inelastic 
scattering 

1 Introduction KTaOs, and strontium titanate, SrTiOa. Both exhibit a 

large dielectric constant increase on cooling connected with 
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The last decade has seen a renewal of interest in the 
paraelectric perovskite crystals [||g] potassium tantalate, 



the softening of a transverse optic zone-center mode (TO- 
mode). The softening of the TO-niode frequency f2sm is 
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both crystals remain paraelectric down to the lowest tem- In addition to these experimental difhculties, the mod- 

peratures (T). These materials, alternatively called incipi- elisation of the phonon spectra in tetragonal SrTiOs is in- 

ent ferroelectrics or quantum paraelectrics, exhibit a wide trinsically more complex due to the presence of the three 

range of anomalous thermal properties associated with the structural soft modes originating from the cubic- phase 

low frequency TO-mode and its mixing at finite wavevec- R25 zone-corner mode. For all these reasons, it appeared 

tor with the transverse acoustic (TA) degrees of freedom preferable to first investigate in detail the ferroelectric TO- 

l^jQ. Of particular interest is the observation of an extra mode and acoustic mode dispersions in KTaOs where the 

Brillouin scattering doublet with a sound-like dispersion quantum paraelectric state is not perturbed by the pres- 

relation, originally attributed to second sound P,pinl- This ence of other low- lying optic branches. 

new excitation appears on top of the broad quasi-elastic t r ^- ^i i r -^ ^- ■ r^rn r\ 

''^^ ^ ^ Iniormation on the low frequency excitations m KiaUa, 

central peak first reported by Lyons and Fleury Ipl. Un- . -i 1 1 r ^^ • • ^ 1 n 

t-i' IS available trom scattering experiments such as, Kaman 

derstanding these low-T features requires a complete and g^ hyper-Raman ^,^, Brillouin §gg, and neutron 

accurate knowledge of the low- frequency phonon disper- , , . iRi Fn;i iToi FT?] Rrl j -c ■ c j * 

° -1 J ^ f scattering p, |17| , |18| ,|19|,pC|, and trom infrared spectroscopy 

sions not only along the principal axes but also in non- rPTTllo^ tj- ^--c i- • 1 -lui-c j-i 

■' '^ t- '^ [^|2^. Indirect information is also available trom dielec- 

■^ ^ ■ trie [g^ and thermal |£J| measurements. Neutron-scattering 

studies have revealed a strong anisotropy in the TA-mode 

mi r o ^TTr^ • i- ^ i ii ii i r dispcrslou, suggcsting an anisotropic TA-TO interaction 

ihe case of SriiOs is more complicated than that of ^ ' *"= ^ ^ 

T^rr r^ u c rr-n. j 4-u u 1 u- mechanism l3,M. Perry et al. PGl measured the complete 

KiaOa, because briiOa undergoes the well-know cubic- u'l — ' •' L— i' ^ 

. ,, 1 ,.r i-ij_- i i li -i- i set of phonon dispersions along the three principal direc- 

to-tetragonal antiterrodistortive structural transition at '^^ ^ ° ft- 

Ta ^ 105 K |,0,|ll[li, below which the unit ceU is dou- *^°'^' ^* '^^^''^1 temperatures between 4 K and 1200 K. 



, , , ^ . i ii • , 1 i 1 • ii 1 • i i 1 1 The data were analyzed using a self-consistent anharmonic 

bled. Owing to this, the crystal m the distorted phase ^ 

rrn ^ m \ r i i rj • ii -i-i shcll-model taking account of the anisotropic non-linear 

[1 < la) can form two types of domains: three tetrago- ° 

,1 . -.1 -1 • n 1 . •;! r- ; 1 polarizabihty of the oxygen O^^ ion fEol|25|,|26|l. 

nal domains with the c-axis parallel to either one of the ^ ■' ■"' iz-il—l'l—if 

three <100> cubic directions and, in each case, two an- In the present paper we are actually interested in a sim- 

tiphase sub-domains corresponding to opposite rotational pier phenomenological approach limited to the low energy 

amplitudes for the oxygen octahedra. While the tetragonal branches near the center of the Brillouin zone (BZ). For 

domains can be aligned by application of uniaxial pressure this purpose an expansion of the phonon eigenfrequencies 

|0,pj, the antiphase domain structure is more difficult to (uj) in powers of the wavevector (q) can be used, and such 



control and there is growing evidence that it may give rise a theory has already been developed by Vaks |E7,28|. We 



to anomalous spectral features |T3]. shall show that the model parameters estimated by Vaks 
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pq do not describe the low-energy low-T phonon disper- stat, with its vertical axis along a [001] or [Oil] direction. 

sion curves in KTaOa with the degree of accuracy needed It should be remarked that the thickness of 0.2 cm is near 

for computational purposes. In order to reliably evaluate optimal for cold neutron work due to the large absorption 

the expansion parameters in the model, more extensive cross-section of Ta ((Ta (natural Ta) ~ 40 bn at 6 meV). 

and more accurate data are in fact necessary. Although 

high resolution data from cold neutron measurements are 2.2 Scattering technique 

available for Li and Nb doped KTaOa [g^,^,31,3^, no 



such data exist for the pure crystal. For this reason addi- 
tional neutron scattering measurements were undertaken 
as described in Section || below. The results of these mea- 
surements are presented in Section 0, while Section H de- 
scribes their parameterization with the Vaks model. The 
quantitative low-T multiphonon computations based on 
the Vaks-model parameterization will be published sepa- 



rately 1 33, 3 



The experiments were performed on the cold neutron three 
axis spectrometer IN14 at the Institut Laue Langevin, 
Grenoble, France. The instrument is equipped with a py- 
rolitic graphite (002) flat monochromator, and the scat- 
tered neutron wavevector kf was held fixed at either 1.785 
A~^ or 2 A~^. CoUimation was achieved with 40' and 60' 
SoUer's, before and after the sample, respectively. As an- 
alyzer, a horizontally focusing silicon (111) crystal was 
used in order to suppress the second order neutrons (2A:/). 
The horizontal analyzer curvature broadens the elastic 
2 Experimental conditions and data analysis q-resolution, but does produce a significant intensity gain, 

with only slight effects on phonon line-widths. The in- 
2.1 Sample 

strumental (q, co) elastic resolution ellipsoid was mapped 

The Ultra High Purity (UHP) crystal used for this exper- out using the KTaOs (200) Bragg reflection. It is well 

iment was grown by one of us at the Oak Ridge National represented by a 4-dimensional gaussian distribution with 

Laboratory using the "modifled spontaneous nucleation" widths (hwhm) Aq^ = 0.017 rlu, Aqy = 0.013 rlu, Aq^ = 

method in a slowly cooled flux |P5|. This technique yields 0.05 rlu, and AE = 0.17 meV, where the indices x, y 

single crystals with blue semiconducting parts surrounded and z refer to longitudinal [100], transverse in-plane [010] 

by colorless UHP-material. The sample used here is a 2.5 and vertical [001] axes respectively. The reciprocal lattice 

X 1 X 0.2 cm^ polished colorless plate, with faces nor- units (rlu's) are defined as ^ = qa/2n, where a — 3.983 A 

mal to the <100> principal axes. The effective mosaic of is the lattice parameter of KTaOa. 

the plate, originally below 5', increased to about 20' in Excitations were measured on the Stokes side (neutron 

the course of the experiments, possibly owing to repeated energy loss) over a typical energy range from up to a 

thermal cycling. The crystal was mounted in a ^He cryo- maximum of 8 meV. The phonon branches along the 3 
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principal directions were scanned up to ^ ~ 0.5 rlu with Jj, 



special attention given to the vicinity of the BZ center, 5,- fa uj) ~ — ?; — ^^ — ^^-^ (2) 



1 [l + n{uj)]ojrj{(i) 



1^1 < 0.2 rlu. A total of approximately 150 high resolution 
inelastic scans was obtained at 10, 24, 39, 60 and 80 K, 
amounting to an equivalent of 4 weeks measurement time. 



The actual fit is obtained using a least square convergence 
criterion implemented in a Marquardt-Levenberg gradient 
algorithm |39| . A typical energy scan along a principal 
axis is illustrated in Fig. H. A typical off-axis wavevector 
scan, perpendicular to a principal direction, is shown in 
The inelastic neutron scattering data treatment was per- ^ig. I As already remarked in |I|,|o|, the asymmetrical 
formed using the usual procedure in the case of anisotropic 1^=^^ P^^^^ ^^^^^ i^ ^ig' ^^ ^^1^^ accounted for by the 
dispersion surfaces HM- First, the gaussian 4- finite wavevector resolution combined with the steep and 



2.3 4D-Monte Carlo adjustment 



dimensional (q,c.) modelization of the resolution elhp- anisotropic dispersion of the TA-mode. 
soid was performed using the Popovici method ||36| im- 
plemented in Rescal [ p7[ . The model ellipsoid was then 3 Experimental results 

convoluted with a locally anisotropic quadratic expansion 

3.1 Zone-center soft TO-mode 

of the dispersion surface using a Monte Carlo algorithm 

around a reference point (^0,^0)- The expansion is taken frenuencv H and damning T of the soft F, 

(r'15) TO-mode at the BZ-center were obtained as a func- 

, /2/c c \ ( , \ „ c \'^ \ ( ^ c ^2 I /„ c ^2 n\ tion of temperature. As shown in Fig. pi the squared TO- 
uj (t-4o) = (^o+c//4//j +(c±.i4±.ij +(ci,24±,2j , (Ij ^ b a -i 

frequency follows a standard Curie- Weiss law above ^ 30 

where ^/ / is the projection of the reduced wavevector shift, 

K, and saturates at a finite value i?^^ ^ 6 meV below 

C — Co: along a direction set by the symmetry of the dis- 

10 K. This is in perfect agreement with previous hyper- 

persion surface at the reference point ^0, while ^^^1 and 



Raman determinations [^|1^,^. For the mode-damping 

^_L 2 are the components of the shift normal to this direc- l-. 

coefficient (Fig. ^, the present neutron-scattering mea- 

tion. The c's are appropriate coefficients defining the local 

surements also agree perfectly with recent hyper-Raman 

curvatures of the dispersion surface, whose numerical val- , , , r—. 

data on a sample of similar quality ||4C|| . 

ues are obtained iteratively. Expression (1) is of course 

modified appropriately if the reference point is at ^0 = 0, 

3.2 Low energy phonon-dispersion curves 

in which case loq is a local minimum for the TO-mode. 

Phonon line-shapes are taken as damped harmonic oscil- Fig. shows the phonon-dispersion branches measured 

lators ||38| with eigenfrequency ujj and damping coefhcient in the high symmetry directions. These data complement 
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the previous results from Refs. p.h9,M. The transverse and Shirane first studied this anisotropy at 20 K and 300 K 

acoustic phonons measured at Q = [^, 0, 2] and by mapping the LA, TA and TO-niode dispersions in the 

[C/\/2, C/\/2, 2], are polarized along the [OOi] axis. They vicinity of the [100] axis. They found that the TA and TO- 

both correspond to the elastic constant C44 at small ^ and mode dispersions form valleys that become deeper as T 

hence have the same dispersion over the central region of decreases, while the LA-mode exhibits a weak anisotropy. 

the BZ. The coupling between the TA and TO-modes low- They ascribed the characteristic diffuse streaks observed 

ers the acoustic phonon frequency as the TO-mode softens by X-ray scattering |4^ to the strong anisotropy of the 

on decreasing T y,0- This feature affects all TA branches TA-mode dispersion. 

polarized along <100> directions. Hence it is not seen Our off-axis scans confirm this behaviour for the TA 
on the TA-branch propagating in the [111] direction, nor and TO-branches around the [100] and [110] directions 
on the TA-branch propagating along [110] and polarized (Figs. ^ and 0). The energy accuracy of ~ 0.1 meV which 
along [110] (not shown in Fig. @). is achieved on-axis deteriorates slightly as one moves away 
The slopes of the TA-dispersion curves in the vicinity from the high symmetry direction. It is clear from Fig. |7| 
of the BZ-center are in good agreement with the elastic that the Vaks model using the original KTaOs parameter 
constants that are known from neutron and Brillouin scat- set [Eq] is not satisfactory for an accurate parameterization 
tering [P,EOlM. In contrast to the case of the TO-mode, of the dispersion surfaces away from the principal axes, 
no measurable damping could be extracted from the TA even at q values as low as 0.05 rlu. 
and LA-mode experimental line-shapes. This sets an up- 
per Hmit for Fta and Fla at about 0.05 meV. Finally, 4 Parameterization of the phonon spectrum 

as seen in Fig. |g, the original Vaks model parameters [g8[ 

Two alternate approaches have been used to describe the 
do not describe satisfactorily low-T TA-modes along high 

phonon spectrum of KTaOs. First, the phenomenologi- 
symmetry axes for |q| > 0.1 rlu (w > 2 meV). 

cal description due to Vaks [^^'Hl] consists in expanding 

the lattice Hamiltonian in powers of q, taking only into 
3.3 Anisotropy of the spectrum 

account the low lying phonon branches, and the interac- 

In order to characterize the phonon-spectrum anisotropy, tions anrong them, as allowed by symmetry. The expan- 

some off-principal axis scans were performed. Starting from sion, limited to the 5 lowest-energy excitations (two TA, 

a reciprocal space position on a high symmetry axis, one LA, and two TO-modes), is actually equivalent to a 

q-scans perpendicular to this axis have been performed. Ginzburg-Landau free energy expansion with respect to 

A typical q-scan at constant oj is shown in Fig. 0. Some dielectric polarization and elastic deformations. Similar 

w-scans at constant off-axis q were also performed. Comes ideas have been discussed by Axe et al. [g| in order to 
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account for the strong T-dependent depression in the TA- nation in order to obtain a rehable description of low- lying 

mode dispersion owing to the repulsion with the softening dispersion surfaces. 

TO-mode. In the Ginzburg-Landau formalism, this effect 

is described as a symmetry-allowed interaction between 4_j Effective Hamiltonian 

the shear strain and the gradient of the electric polariza- 

One considers first the harmonic part of the lattice Hamil- 
tion. ^ 

tonian 7i of a cubic perovskite crystal. Let Ui(q) be the 

A more ambitious approach was developed by Migoni spatial Fourier transform of the displacement from equi- 

et al. [g5|. This is a microscopic model which aims at librium for the ith atom in the unit ceU. The mass of this 

describing the complete phonon spectrum across the entire atom is denoted m^ and its effective charge, e^. One finds 

BZ, and which accounts for the softening of the TO-mode [p8l|43,E3|] 
through the non-linear polarizability of the O^^ ion. Perry , 

et al. [gO[ used this model to describe the phonon branches i.q 

of KTaOs along the high symmetry directions. Their ex- + o ^ Un(q)^n,'2(q)uJ2(-q) 

il,i2,q 

tensive neutron measurements allowed to determine the V^ /_at7i/ _.\ , ^ 

full set of 15 model parameters. However, the resolution 
of the experiments, and thus the accuracy of the resulting 
model, are not better than ^^0.5 meV. 



2 

»i,»2,q 

^e.Ui(q)E(-q) + ^^|i^(q)|^ (3) 

i,q 



where eoo is the high frequency dielectric permittivity of 

the crystal. The tensor <P accounts for the "short-range" 

part of the potential energy, i.e. for the contribution of the 

At low temperatures, the thermodynamic and phonon- short-range forces plus the analytical (in q, at q — > 0) part 

kinetic properties of KTaOs are controlled by the 5 lowest of the dipole-dipole interaction. The non-analytical part is 

modes included in the Vaks model [ p7| , [2q ]. Since our aim taken into account by introducing the Fourier transform 

is not to explain the origin of the soft mode, but rather of the macroscopic electric field E(q), which is given by 

to use the phonon spectrum to account for the low tem- Poisson's equation, 

perature anomalies observed in the quantum paraelectric 

q • [eooE(q) + P(q)] = . (4) 

state, the Vaks parameterization is a very convenient one. 

It includes 8 adjustable parameters, but 4 of these are The Fourier transform of the macroscopic electric polar- 
well determined from light-scattering data. As the physical ization is defined by 
basis of the model is only detailed in Rcf. [pSJ, we sum- -, 
marize it below, specializing to the cubic case which is of '^ i 
interest here. We shall then refine the parameter determi- where Vc — a^ is the unit-cell volume. 
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Using a routine procedure, Eqs. (0), (0), and (h) give from the long-range part of the dipole-dipole interaction, 
an effective Hamiltonian that can be diagonalized in terms ,^e \ V^ f zi,^aa'Zh\i3i3'\ ,„\ 

(^b,6')a/3= 2^ ^3 Ina'nfB', (9) 

, „, \ £00" / 
of normal phonon coordinates. However, following Vaks 

where n = q/|q|, and Zb^aa' = J^i^i^iaM' /V^ is the 
p7|,p8| the most convenient form of the Hamiltonian for 

effective charge tensor for the 6-set of optic modes in the 
the present purpose is obtained using a linear transforma- 

limit q — > 0, which in view of the cubic symmetry is oc Saa' 
tion of the displacements Ui(q) which diagonalizes only 

the short-range part of H in the limit q -^ 0. One selects 

One is only interested in the low frequency part of the 

a particular u to describe the acoustic modes by 

spectrum. Following Vaks P^Eq], one first diagonalizes 

ui(q) ~ u(q) . (6) the optical part of (|h in the limit q ^ 0, keeping only 

the terms containing the restoring short-range force con- 
One defines then the optical displacements Xf,(q) by 

stants of the optic modes , Ab, and the non-analytical term 

1 ^ 
Ui,a{<l) — ui,a(q) — } Lia,bi3Xb,p{o^) , (7) ^f (,' • C'ne fiuds that the latter lifts the degeneracy of the 



mi 

I3,b=l 



optical modes, splitting each set into a pair of degenerate 
where a and P are cartesian indices. In writing (R), one 

TO-modes, whose frequency squared is equal to Afc, and a 
takes into account that the short-range part of the lattice 

longitudinal (LO) mode. One pair of TO-modes in the set 
Hamiltonian of a cubic perovskite at small q contains 4 

b corresponds to the soft mode, which has a frequency f2sm 
sets of 3-fold-degenerate optic modes that correspond to 

much lower than all other optic modes. It is associated to 
the same vector irreducible representation. The subscript 

the hardest LO mode [H, labelled LO4 in the Perry et 
b indexes these sets. With the new variables, Eqs. ^, (^, 

al. denomination ^^. We denote the corresponding Xt, 
and (g) give: 

by A and the corresponding optical displacement vector 

H = - 2_^[u(— q)u(q) + u(— q)A(q)u(q) Xb by x. The anharmonic interactions, not included in a, 

q 

, v^ ( • f \ • I- \ , \ I \ I \ are taken into account with a standard quasi-harmonic 

+ 2^{x6(-q)xb(q) + AbXb(-q)xf,(q) ^ 

approximation. The frequency of the soft mode becomes 

+ 2u(-q)n(q)xb(q)} 

then temperature dependent through the T-dependence of 

+ J2 M-Ci){Sb,b' (q) + i>lb' (q)}xb' (q)] , (8) 

b^b' A = f^sm- 1^ futher considerations, only the subset consist- 
where Af, are the restoring short-range force constants of ing of the two soft TO-modes and of the 3 acoustic modes 
the optic modes. The tensors A, Vb, and Sb^b' correspond is kept. One neglects the other optical modes as well as 
to the contributions of the short-range interactions. These their interactions with this subset. As shown by Vaks [ p7| 
analytically go to zero as q^ for q -^ 0. The tensor S^^^, pa, the accuracy of this approximation for the description 
is non-analytic at q — > 0. It represents the contribution of a phonon of frequency to is controlled by the small pa- 
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rameter w^/J7^ , where fiop is a typical hard optical mode The components of the tensors A, S and V that give non- 
frequency such as the hard LO4 mode. Hence for the part vanishing contributions in the subspace (cci, 2:2, ui, ^2, U3) 
of the spectrum that we are interested in (w ~ f^sm), it in the small-g limit are 
results from the Lyddane-Sachs- Teller relation that this 



parameter is of order eoo/eo, where to is the static dielec- 
tric permittivity of the crystal. 

Technically, the separation of the Hamiltonian contain- 
ing only the 5 modes of interest is done to within the 
aforementioned accuracy by considering the problem in a 
rotated reference frame X' = {X[, X2,X^), rather than in 
the usual cubic frame X = {Xi,X2,Xs). The axis X'^ is 
taken parallel to q. Specifically, the rotated frame can be 
defined by the relation X' = il/X, where the matrix M is 
given as a function of the direction of q (n = q/q) by 



A^q^iAig' + Atg^+Aag'^) 
S = q\Stg' + Sag""), 
V^q^Vtg' + Vag''), 



(12) 



with 



J 

?Q/3 



lani3, 



JaP = ^ap - rianp, 



Jaf3 = laPfSn.yns. 



(13) 



Af(n) 



n± Til 

y n2/n_L -niU'i/n^n'i I 



(10) 



The tensor jap-yS is defined, in the cubic frame, by jap-yS = 
liia — f3 = "/ = S, and jafs-yS = otherwise. In that 
frame. A, S, and V stand for acoustic, optic, and acoustic- 



where n is written in the cubic reference frame X, and ^ptic coupling tensors respectively, split into transverse g\ 
nj_ = I — Uj^. longitudinal g', and anisotropic g"^ parts, as seen in ([l2|). 

One sees that in the rotated frame the components The coefficients Si and Vi do not appear in (|l|) since x is 
of the acoustic and optical bare mode displacements (as taken as purely transverse, 
defined in (g) and (0)), associated with the 5 modes of in- 
terest are just xi,X2jUi,U2, and U3, while 3:3 is essentially 
the displacement of the LO4 mode which has been elim- 
inated. Finally, the 5-mode Hamiltonian can be written 



Since the kinetic energy terms are already diagonal, 
the phonon-dispersion curves for the 5 branches are given 
by the square roots of the eigenvalues of the following 
matrix: 



q 

-I- AX_gXq + x_gS'(q)xg (11) 

+ 2u_5t>(q)x,] . 



■H'-'H<D^'H,si\cil\,St,At,VuAi) 

+ Hamsiq,Sa,Aa,Va). (14) 
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The isotropic matrix, Tiis, is 
/ 



V 



\ + Stq^ Vtq^ 

\ + Stq^ Vtq^ 

Vtq^ Atq^ 

Vtq^ Atq^ 

Aiq^ I 

and the anisotropic one, "Hanis, is 

' Sahii Sahi2 Vahii Vahi2 K/il3 ^ 

Sahl2 Sah22 Vahi2 14/l22 K/l23 

VJlii Vahi2 Aahii Aahi2 Aahis 

Vahi2 Vah22 Aahi2 Aah22 ^q^23 

\Vahi3, Vah23 ^a^llS ^a^23 Aah^s. J 

in which 

n 2 2 

hii = -4^ , 



(15) 



(16) 



tis 



_ n2ri3 /'„2 



^K-«^) 



hi: 



_ niri2ri3 / 2 



^(r.^ - ni) , h2, - ^(n^ni - n| - nf) 



/i22 = 2n^(n^ - ^^) , /i33 ^n\ + n\ + n\ . 

The eigenvalues (^i(q) of 7i(^)(q) are the squared phonon 
frequencies, and the diagonalization matrix Lij(q) reflects 
the mode mixing. Thus, the renormalized acoustic u and 
optical X displacements in the canonical cubic axis system 
are obtained through a linear combination of bare mode 
displacements: 



^ M,,4fq' and u,. q = Y, M,v_2l2 



(/) 



J^„t = l,2 



iac=3,4,5 



(17) 



where Mj indicates the column j in matrix M . 

In the matrix Tiis, the parameters St, At, and Ai ap- 
pear as squared dispersion-curve slopes whereas Vt is a 
TA-TO coupling parameter, and 5*0, Aa, and Va contribute 
to the off-principal axis squared slopes. 



4.2 Adjustment of the phonon-dispersion surfaces 

We now use the above model to parameterize the mea- 
sured phonon spectra. Since the parameter set given by 
Vaks for KTaOs (see Table |] and Figures |^ and ^ leads to 
typical deviations of 0.5 meV or more with respect to our 
inelastic neutron scattering measurements, it is necessary 
to recompute the parameter values. The model contains 
8 parameters: St, At, Vt, Sa, Aa, K, Ai, and A = i?^^, 
where only the latter is appreciably T-dependent. Three 
parameters can be expressed in terms of the elastic con- 
stants Cij and the density p = 6.967 g/cm^ of the crystal 

Al = (Ci2 + 2Cii)/p , 

At = Cii/p , (18) 

Aa = [Cn - Ci2 - 2Cii)/p . 

The values of the elastic constants are known from Brillouin- 
scattering experiments ||,|4|: Cn = 4.344 ±0.049, C12 = 
1.027 ± 0.038, and C44 = 0.990 ± 0.020 in 10" N/m^. 
The T-dependent values of A are available from the hyper- 
Raman measurements 115100, Eo] . The values used presently 
are given in Table 0. Thus, there remain only 4 adjustable 
parameters at our disposal, namely St, Vt, Sa, and Va- 
We adjusted our data on the 5 lowest frequency modes for 
the <100>, <110>, and <111> directions to the square 
roots of the eigenvalues of the matrix (nj) . The fit was per- 
formed with a Marquadt-Levenberg algorithm and a least 
square criterion, for the temperatures, 10, 24, 39, 60, and 
80 K, keeping the parameters temperature independent. 
The overall correlation coefficient when using 151 neutron 



10 E. Farhi et al.: Low energy phonon spectrum and its parameterization in pure KTaOa below 80 K 

and Brillouin measurements is r^ = 0.995, and the deter- resolution ellipsoid was used in order to extract phonon 
mined model parameters are given in Table |l|. This table frequencies and dampings from the experimental spec- 
also contains the sets of parameters determined by pre- tra. The present results improve significantly on previ- 
vious workers |^,^,^. Appreciable difi^erences in some ous phonon studies in KTaOs at low-T [p,17;lq,p^,p0[. 
of these can be seen. The obtained parameterization of In particular, the accuracy on the phonon frequencies if 
the spectrum has been tested by comparing the phonon of the order of 0.1 meV along symmetry directions and it 
frequencies measured off principal directions with those appears to remain better than 0.3 meV in non-symmetry 
calculated from the parameters, as shown in Figs.0 and ||. directions. The phenomenological model by Vaks P7| , P8| , 
One sees that the model gives also a reasonable descrip- which contains 8 parameters (4 of them being determined 
tion for the off-axis dispersion. from light-scattering experiments) was used for the pa- 
The model provides a global description of the low- rameterization of the low-energy phonon spectrum. By 
frequency spectrum of KTaOa as illustrated in Fig. || and adjusting 4 temperature-independent parameters an ac- 
|lO| . A remarkable feature of the phonon spectrum is the curate description of the spectrum was achieved over the 
strong anisotropy of the TO and TA modes. We conclude temperature range from 10 to 100 K. 
that the The relatively simple description of the temperature- 
Vaks model provides a very reasonable parameterization dependent spectrum of KTaOs which is thus obtained pro- 
of the five lowest phonon modes of KTaOa {ui < 8 meV) vides a convenient basis for quantitative calculations of 
in the central part of the Brillouin zone (|q| < 0.25 rlu) all the properties that are controlled by the low frequency 
and in the low temperature region (10 < T < 100 K). phonons. In particular it becomes possible to calculate 
The achieved accuracy is '^ 0.1 meV in the vicinity of the thermodynamic and transport properties at low temper- 
high symmetry axes, and remains better than ~ 0.2-0.3 atures, such as phonon-density-fluctuation processes and 
meV away from these i.e. of the order of the experimental second sound, in the low-T quantum paraelectric regime, 
uncertainties. Such calculations will be the object of a forthcoming pub- 
lication. 
5 Conclusion The authors thank Alain Brochier for his efficient tech- 



We presented the results of cold neutron high resolution 
inelastic scattering measurements of KTaOs at low tem- 
peratures, for the lowest phonon-dispersion branches and 
for wave vectors |q| < 0.25 rlu. A Monte-Carlo convolution 
technique based on a modelisation of the 4-dimensional 



nical support during the IN14 experiments at ILL. 
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r(K) 


80 60 39 24 10 


A {meVf 
(uncertainty) 


30.8 23.1 14.9 9.2 6.2 
(0.5) (0.8) (0.5) (0.4) (0.3) 



Table 2. Square of the soft mode frequency. A, used in our 
refined Vaks parameterization. 
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Iq^^l = 0.05 rlu along [110] 
T = (002) 
-ffl =1.27meV 
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Energy (meV) 



Fig. 1. Lj-scan : TA phonon in the [110] direction (|q| = 
0.05 rlu), polarized along [001]; KTaOs, T = 15 K. The fit is 
performed using a '4D' Monte Carlo convolution technique. 
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|q^^| = 0.04 rlu along [110] 




1.96 1.98 2 2.02 

Qj^(rlu) along <001> 



2.04 



Fig. 2. q^-scan across a TA 'valley' around |q] = 0.04 rlu 
in the [110] direction; KTaOs, T = 10 K. The fit is performed 
using a '4D' Monte Carlo convolution technique. 
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Fig. 3. KTaOs: soft TO-mode squared frequency at the zone 
center compared to previous hyper-Raman measurements mt . 
The mean energy uncertainty for the neutron data is AE ~ 
0.1 meV. 
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Fig. 4. KTaOs: Damping coefficient (half width) of soft zone 
center TO-mode compared to previous hyper-Raman measure- 
ments |40|. The mean uncertainty for the neutron data is AF ~ 
0.1 meV. 
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Fig. 6. KTaOa; Dispersion curves measured at 10 and 39 K 
in the [100] direction; the original Vaks parameterization pq] 
is shown as dashed lines. The modified Vaks model (this work) 
is shown as solid lines. 
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KTaO3:|^| = 0.05rlu,T=10K 



30 20 10 10 20 30 

(deg) at |^| = 0.05 rlu 



Fig. 7. Measured dispersion curves in KTaOs at 10 and 39 
K near the [010] direction in the (001) scattering plane: Q — 
[(2+5 sin S), (^cosS), 0]; The dashed and solid lines refer to the 
original [Eq] and modified (this work) Vaks parameters. 
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Fig. 8. Anisotropy of the TO phonon branch near the [010] 
symmetry axis in KTaOs: Q = [(2+5sin6'), (Ccos6'), 0]. The 
lines are calculated with the modified Vaks parameters. 
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Fig. 9. KTaOs: The low frequency phonon spectrum at T = 
10 K, for q = [(^ sin 61), {£, cos 6 / y/2) , (£, cos 6 / y/2)] with |5| = 
0.05 rlu, as computed by our Vaks-type parameterization. The 
directions [110], [111] and [001] correspond to the angles 0° or 
180°, 35.3° or 144.8°, and 90°, respectively. The experimental 
data points are shown as open circle. 
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Fig. 10. Highly anisotropic coupled acoustic and optical 
phonon surfaces at T = 10 K in KTaOs showing a deep valley 
around the [100] direction; calculated from the modified Vaks 
parameters. 
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{meV/rluf 


St 


At 


Vt 


Sa 


Aa 


Va 


Al 


Ref. HH 


4860 


1836 


2808 


10800 


2160 


-5724 


4320 


Ref. § 


4670 


1830 


2800 










This work 


4828 


1553 


2451 


8956 


2265 


-3087 


4552 


(uncertainty) 


(20) 


(10) 


(30) 


(400) 


(40) 


(290) 


(30) 



Table 1. Low temperature Vaks parameters of KTaOa, determined from a global fit of neutron and Brillouin data. See Sec. 
4.1 for the meaning of the parameters. 
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Fig. 5. Dispersion curves of KTaOs low energy phonons between 10 and 80 K along the high symmetry axes. The mean energy 
uncertainty is 0.1 meV. The LA mode is measured at Q = [0, 0, 2+^]. The lines correspond to calculations using the Vaks 
parameterization presented in this work. 



